TIFR/TH/10-01 



O 

(N 
O 

Q 

(N 



I 

or 
43 . 



> 
On 

cn 

o 
o 



Lattice QCD predictions for shapes of event distributions 
along the freezeout curve in heavy-ion collisions 

R. V. GavaQ and Sourendu GuptcQ 

Department of Theoretical Physics, 
Tata Institute of Fundamental Research, 
Homi Bhabha Road, Mumbai 400005, India. 

We present lattice QCD results along the freezeout curve of heavy-ion collisions. The variance, 
skew and kurtosis of the event distribution of baryon number are studied through Pade resumma- 
tions. We predict smooth behaviour of three ratios of these quantities at current RHIC and future 
LHC energies. Deviations from this at lower energies signal the presence of a nearby critical point. 

PACS numbers: 12.38.Gc, 25.75.-q, 11.15.Ha, 05.70.Fh 

All matter that we know undergoes phase transitions as external conditions change. Strongly interacting matter, 
described by QCD, is not thought to be exceptional. Although the phase diagram of QCD matter is not yet established, 
it has been the subject of intense theoretical and experimental scrutiny in recent years. Finding a critical point in 
this system would be major landmark. One idea is to look at quantities which are non-monotonic near a critical 
point. In this letter we give lattice predictions of a few such observables; we give precise results for the backgrounds. 
These are the first lattice QCD predictions of any quantity along the freezeout curve of heavy-ion collisions. Any 
non-monotonic behaviour over these is a signal for the critical point; and the lattice results near the critical point do 
show such behaviour. 

Twenty years ago it became clear that QCD has no finite temperature (T) phase transition since all quarks 
are massive, as evidenced by a non-vanishing pion mass. Lattice computations showed that there remained a fairly 
abrupt change in thermodynamic properties which could possibly influence the physics observed at the RHIC and 
LHC. About a decade back it was realized that QCD could nevertheless have a critical point at finite baryon chemical 
potential (hb) [21 ■ There is evidence for such a critical point from modern lattice computations with two light 
dynamical flavours [H, 0] , and from an earlier small volume computation [f| . The behaviour of QCD with two light 
and a single heavier flavour is not yet settled, although evidence for a critical point is mounting [6(. 

One question of interest is whether this critical point can be found in experiment. The results of [H, 0] suggest that 
heavy-ion collision experiments at moderate colliding energies, achievable at RHIC, FAIR and NICA could find it. In 
these collisions of heavy-ions, a fireball is created which evolves chemically and thermally before freezing out. The 
primary question of experimental interest is how close the freezeout is to the critical point. In this work we employ 
Pade resummation of lattice results to obtain predictions which can be directly compared with experiments. 
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FIG. 1: The critical point of QCD is the end point of a line of coexistence of two phases. One of the interesting questions for 
heavy-ion physics is how close the freezeout curve is to the critical point, i.e., which of the three is the true position of the 
freezeout curve. Our results, discussed below, show that T c — 170 corresponds to a line like A, whereas an increase of 5 to 10 
MeV in T c could shift the results to the curve C. 



There has been some discussion on suitable observables for the identification of a critical point. The first suggestion 
was to study the width of the distribution of event-by-event (E/E) observations of momenta of particles [2j. Later 
it was realized that a measurement which was suggested in a different context could be suitable: that of the width 
of E/E distributions of conserved quantities Q- For the baryon number, this is proportional to the baryon number 
susceptibility (BNS) [1, and is expected to diverge in a thermodynamically large sample in equilibrium at the 
critical point. It was soon pointed out that non-linear susceptibilities (NLS) have stronger diver genc e As 
a result, higher cumulants of the E/E distributions are als o g ood signals for the critical point llj, [12J. It has been 



pointed out recently that ratios of cumulants are better [13j . since ill-determined parameters such as the fireball 
volume do not appear, and they are directly comparable to predictions from lattice QCD — 

T X {3) (T,n B ) T*x W (T,»b) T X M(T,Hb) m 

X (2) (T,hb) x K) {T,ijlb) x x >\T,11.b) 

where X (T,fiB) is the n-th order NLS, obtained by taking the n-th derivative of the pressure with respect to fiB- 
As one can see, the ratios are not all independent, since m,2 — mim^. These ratios are also of interest in lattice 
studies of the phase diagram of QCD, since they provide estimates of the distance to the nearest singularity of the 
free energy. 

At the critical point, (T E , fig), the BNS, x^(T E , Ms): diverges. Since this critical point is in the Ising universality 
class, the eigendirections of the renormalization group scaling transformations can be classified as a "thermal" direction 
t and a "magnetic" direction h, as illustrated in FigureQ] The baryon-baryon correlation length diverges with exponent 
y t as one approaches the critical point in the direction t, and with exponent yh along h. Once these exponents are 
known, the divergence of the BNS can be worked out. In the present work we do not need the value of the exponent 
8, which controls the divergence of the baryon number susceptibility through x (/Ab) — Ims — Msl* 5 - All we need to 
know is that 

X^l^-^ir"- 2 . (2) 

As a result, the ratios of successive NLS diverge as 1/|/zb — Ms I as one approaches the critical end point. Due to the 
CP symmetry of QCD, a critical point at Hb implies another at with the same critical exponents. Hence, the 
divergence of the ratios of successive NLS can be rewritten as l/|/i 2 — (mb) 2 |- 
In lattice calculations one uses the Maclaurin series expansion of the BNS, 

= + ^X {A \T)z 2 + ^TV 6 >(T)z 4 + irV 8 )(T)z 6 + 0(z% (3) 

organized as an expansion in the dimensionless quantity z = (J.b/T at fixed T. The Taylor (Maclaurin) coefficients, 
T n - 4 x (n) {T), which are also dimensionless, are evaluated through lattice computations at /is = 0. We have written 
the expansion out to the order up to which the coefficients are currently known [3, 0] ■ The successive estimators of 
the radius of convergence of this series are 



\ (n) (T) 

<n +2 = (n + 2)(n + l) T2 ^ 2 / (T) . (4) 



These have been used to estimate the position of the critical point in 0, 0] , where the finite volume effects in these 
ratios were investigated. It was found that r 2 „ +2 for n < 6 were insensitive to the volume, provided V > (4/T) 3 . 
In terms of these quantities the expansion of eq. ([3]) can be written as 



J>2 J^2 



1 + — + 2 2 + 2 2 2 + 0(z 8 ) 

r 24 r 24 r 46 r 24 r 46 r 68 



(5) 



Clearly, the ratios of eq. (fT]) depend only on the three quantities r| 4 , r| 6 and rg 8 , and not on x (T)- By taking 
formal derivatives of the expansion above one obtains the higher order susceptibilities. Using these one arrives at the 
formal series expansions of the measurements we are interested in — 
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It has been demonstrated that series expansions are not a good way to extrapolate physical quantities to finite 
chemical potential, and better techniques, such as Pade resummations are needed Since we know that mi and TO3 
have simple poles and 1TI3 has a double pole, Pade approximants of appropriate order are particularly suited to this 
problem- 
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zP^z 2 ; a, b), m 3 = -Pf(z 2 ;a', b') 



m 2 



p 2 \z 2 - 



a, b), 



(7) 



where the notation P^ indicates that the numerator is a polynomial of order L and the denominator is of order M 
and P%m indicates that the denominator is the square of a polynomial of order M. Since we have three terms of each 
series at our disposal, in each case we can use L = and 1. When longer series expansions become available, one can 
use larger values of L. The parameters clq ■ ■ ■ <Xl-i, af - ■ ■ a' L _ 1 , and anj - • 5l-i in the numerator and b\, b' x and 61 in 



the denominator are obtained by matching to the series expansions [14j. This is done within jack-knife blocks. We 
check that the relation m 2 — minis is satisfied within errors. 



V Snn 


z 


T/T c 


P(N t = 6) 


P(N t = 4) 


5.0 


4.7 (2) 


0.70 (3) 




5.20 


7.7 


3.02 (1) 


0.82 (2) 




5.22 


11.5 


2.08 (7) 


0.89 (1) 


5.39 


5.24? 


18.0 


1.39 (5) 


0.94 (1) 


5.41 


5.275 


19.6 


1.29 (4) 


0.94 (1) 


5.41 


5.275 


27.0 


0.96 (2) 


0.96 (1) 


5.415 


5.275 


39.0 


0.68 (2) 


0.97 (1) 


5.415 


5.275 


62.4 


0.44 (2) 


0.97 (1) 


5.415 


5.275 


200.0 


0.142 (5) 


0.98 (1) 


5.415 


5.275 


850.0 


0.034 (1) 


0.98 (1) 


5.415 


5.275 


2500.0 


0.0115 (4) 


0.98 (1) 


5.415 


5.275 


5000.0 


0.0058 (2) 


0.98 (1) 


5.415 


5.275 



TABLE I: The values of z = hb/T and T/T c at freezeout for various \/S NN are shown for T c = 170 MeV. The numbers in 
brackets are errors in the least significant digit and come from the statistical error quoted in [TJ[. There are also statistical errors 
on the scale setting on the lattice [3, 4]. The combination of these errors may result in the same coupling, /3, corresponding to 
several different central values of T/T c . Missing entries for /3 mean that the simulations in [3,0] do not correspond to these 
conditions) . 



In order to compare lattice predictions with experiment one needs the values of T and [Ib along the freezeout curve 
as a function of the center of mass energy of the colliding nuclei, \J S NN . These were extracted from data in [ToT ]. 
For a comparison with lattice computations it is necessary to convert to the dimensionless variables z = Ub/T and 
T/T c . We have used a crossover temperature T c = 170 MeV consistent with the current best estimates 18|, [l9(. The 
correspondences between lattice parameters and the freezeout points are given in Table U where we have tabulated 
parameters corresponding to the RHIC low-energy run as well as the LHC heavy-ion runs at both full energy and the 
expected energy at the end of the current year. 

Once all non-thermal sources of fluctuations are eliminated from data, the ratios of eq. ([T]) are directly comparable 
to lattice QCD predictions. Lattice predictions have long been used to interpret experimental data from heavy-ion 
collisions; examples being the qualitative agreement of the equation of state of hot QCD matter with experimentally 
extracted values and the agreement of strangeness production parameters extracted from lattice measurements and 
fits to hadron yields. These have been indirect agreements, since a layer of analysis stands between the data and the 
prediction. The main results of this paper allow direct comparison of lattice predictions and data for these ratios 
along the freezeout curve. Agreement of data with these values are direct tests of lattice QCD, and could potentially 
be the first quantitative tests of lattice QCD in a thermal environment. However it is interesting to point out other 
aspects of this argument. 

If the critical point is far from the freezeout curve over a certain range of energy, then mi decreases with increasing 
y/S NN (since z decreases) and m.3 increases. Using these two measurements and comparing with lattice predictions, 
it is possible to estimate the freezeout conditions: T/T c and hb/T. This method is independent of the usual one 
in which hadron yields are interpreted through a resonance gas picture [i~5| . Comparison of the two methods then 
allows us to estimate T c by inverting the argument of the previous paragraph. Mutual agreement of the values of T c 
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so derived at different \/S NN would constitute the first firm experimental proof of thermalization. If this proof holds 
then one also obtains the simplest and most direct measurement of T c found till now. Since such a thermometric 
measurement can be made reliably with data at large V S NN , where is small, it would remain a valid measurement 
whether or not a critical point is found in the low energy scan at RHIC. 

Our numerical results are based on the measurements of the Taylor coefficients in the expansion of cq. ([3]) which 
were reported in d, 0|- Since the data were taken at two different series of lattice spacings, with N t = 4 and 6, one 
can examine the approach to the continuum limit. In both sets of simulations the pion mass was close to physical, 
being about 230 MeV. The spatial volume varied from (2/T) 3 to (4/T) 3 at all lattice spacings, and up to (6/T) 3 on 
the coarser lattice. The present analysis was performed using a jack-knife estimator of the mean and its error with 
5-10 jack-knife blocks. We used only statistically independent measurements, of which we have more than 50 at each 
coupling and volume. 

A test of the method is to use it at the previously determined value of T E /T c and check whether the extrapolations 
give radii of convergence which are compatible with that previous estimate. The Taylor expansion of x around 
z o = V° B I T is 



X (2) (T,z) _ X {2) (T,z ) x (3) (T,z Q ) _ xW(J> o), 
rp2 ~ + T 2! 



(z-z ) + A \' "> (z-z Q ) 2 +O((z~z ) 3 ). (8) 



This implies that 1/mi and 2/ '1113 are estimates of the radii of convergence of the series. When they are computed 
at T E then zq + 1/mi and zq + 2/7773 have to be compatible with previous estimates of [i E /T E , obtained from the 
Maclaurin expansion of eq. <j5j> - We have verified this for both Nt — 4 and 6. 

This argument also allows us to understand systematics due to lattice volume effects and finite lattice spacing, and, 
therefore, the possible effects of improved actions. Away from the vicinity of the critical point finite volume effects 
are expected to be small as long as VT 3 > 4 3 . Finite volume effects are expected to be larger near the critical point: 
indeed this is one way to check whether lattice extrapolations show genuine critical effects. Lattice spacing effects are 
also easy to understand. Since the Maclaurin series for 27713 starts with unity, lattice spacing correction are expected 
to be negligible and small. Indeed, this is what is observed. On the other hand, since the Maclaurin series for mi/z 
and 77i2 start with the terms X (T) / x^ (T) , one expects that a large part of these lattice effects can be subsumed 
into the lattice spacing effects on these quantities. We can correct for the lattice spacing changes between N t = 4 and 
6 using data from [3, 4]. However, there is not yet enough data to extrapolate to the continuum. Note however, that 
777,1 1 z and 777-2 have a common normalization. Hence we normalize to 7772 = 1 at \J S NN — 5 TeV. 

With these normalizations, we find that at freezeout in the highest RHIC energy of \JS NN — 200 GeV 

mi = 0.140 ± 0.008, 777 2 = 0.95 ± 0.06, "73 = 6.99 ± 0.06. (9) 

The errors shown are purely statistical. The remaining systematic error comes from the imprecision in our current 
knowledge of T c . We can estimate the magnitude of this systematic error by changing T c by 5 MeV. For example, 
using T c = 175 we find, using the corresponding normalization, that 

77ll = 0.14 ± 0.01, m 2 = 1.04 ±0.07, 777 3 = 7.3 ±0.1. (10) 

Hence this source of uncertainty changes results by less than 10% at the highest RHIC energy. Closer to the critical 
region the changes can be much larger. For instance, whether 7712 and 7773 are negative at some \/S NN , as seen in 
Figure [2] depends sensitively on the choice of T c — not being visible for the higher value of T c . 

The dependence of mi , 7772 and 7773 on the energy, yS™, along the freezeout curve is shown in Figure [2] Note that 
the lattice spacing dependence is under good control for 1/ S NN > 50 GeV. In this region of energies z = hb/T is small, 
as a result of which the Pade approximants are almost constant and one has mi ~ z, 7772 — I/2 and 7773 ~ z°. Since 
the freezeout temperature here is smaller than T c , one should compare our results with the predictions of resonance 



gas models. These give 7772 = 1 [16l. 117]. which is close to our results. In our normalization an ideal quark gas yields 
7772 = 2/(37r 2 ), which is very different. 

In the region of energy just below y/S NN — 50 GeV the lattice spacing dependence of the results is not under 
good control yet. This is related to the fact that the critical point observed on the lattice 0, H} lies in this region 
and shifts with the lattice spacings used. Due to this shift in the critical point with lattice spacing, the two lattice 
spacings give different predictions in this region. As mentioned above, the choice of T c is also important. In order to 
show the expected behaviour near the critical point, we have also shown the results for a single lattice spacing — the 
smallest currently available. Note that the behaviour of mi is fairly smooth, although a departure from the behaviour 
mi ~ z is clear. Both 7772 and 7773 show signs of non-monotonic behaviour with the Kurtosis becoming negative. Such 
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FIG. 2: The ratios mi, mi and 7713 along the freeze out curve as a function of \/Snn- The panels on the left compare the 
lattice results for N t = 4 and 6. Those on the right show only N t = 6 along with the hadron resonance gas results marked 
HRG, and a simple power law extrapolation of high-energy data, in order to make clearer the non-monotonicity near the critical 
region. 



leptokurtic behaviour is impossible in the hadron resonance gas. It develops in QCD only when the lattice spacing is 
small enough. 

At very small \/S NN one again leaves the vicinity of the critical point, and the lattice spacing effects are again under 
reasonable control. This happens although z — ^b/T is very much larger than unity, and the lattice measurements, 
made at z = 0, have to be extrapolated very far. The reliability of the lattice prediction at large z is due to the 
stability of the Pade resummations. Note that in this region there is a return to the expected behaviour mi ~ z, 
m,2 — l/z and ~ z . 

Since thermodynamic stability requires x' 2 ' to be positive, the signs of m\ and TO2 immediately tell us the signs of 
both and x' 4 ^ Our observations then show that is positive along the freezeout curve, in agreement with the 
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qualitative argument of [T2J . This implies that those parts of the freezeout curve which current lattice computations 
can access lie entirely below the line of first-order transitions. It is a moot question whether the freezeout line crosses 
the first-order line (i.e., the coexistence line of Figured]) at lower energies or whether future lattice computations at 
smaller lattice spacings will show different behaviour. 

The most pleasant aspect of the current results is that the critical point is close enough to the freezeout curve for its 
effect to be visible as non-monotonic behaviour of the ratios. This is clearest for mi because of the relatively smaller 
error bars. However, the effect is also present in 777,2 and 7773. There is currently some uncertainty in the position 
and size of the peak due to uncertainties in lattice determinations of T c . However, all the available lattice data show 
that the critical point lies within reach of the proposed energy scan at the RHIC. We emphasize that the existence of 
non-monotonicity is a general expectation if there is a nearby critical point. The baseline over which the behaviour 
stands is a precise result from lattice QCD which is new. 

In this letter we have examined lattice QCD predictions for the ratios in eq. (JXJ) . On the lattice they are related 
to the radius of convergence of the series in eq. (jHJ, and they are also measurable in experiments [l3j]. Using Pade 
approximants to resum the series expansion of these ratios (which reproduce previous results on the location of the 
critical point) we predict the ratios along the freezeout curve (see Figure [2]). Current lattice results already show the 
following important qualitative and quantitative trends — 



1. Generally mi decreases with increasing y/S NN and 777,3 increases, whereas 7772 has a flatter dependence on 
energy. This general trend is interrupted in the neighbourhood of a critical point where there are rapid and 
non-monotonic changes in all three observables. 

2. At the highest energies (RHIC and LHC) sufficiently accurate measurements of mi and 1773 may be used to 
deduce the value of the crossover temperature, T c . 

3. The skew and kurtosis are both expected to be positive in the range of energies away from the critical point, 
with a possible change in the sign of the kurtosis in the critical region. 

These computations were performed on the Cray XI of the ILGTI at TIFR. RVG would like to acknowledge the 
hospitality of the Department of Physics at the University of Bielefeld. SG would like to thank Bedangadas Mohanty 
and Nu Xu for discussions. We would like to thank Sandeep Chatterjee for tables of mi, 7772 and 7773 along the 
freezeout curve in the hadron resonance gas model. 
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